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ABSTRACT 



Aims. We investigate the behavior of the frequency-centered light curves expected within the standard model of Gamma 
Ray Bursts allowing the maximum electron energy (7max) to be a free parameter permitted to take low values. 
Methods. We solve the spatially averaged kinetic equations which describe the simultaneous evolution of particles and 
photons, obtaining the multi-wavelength spectra as a function of time. From these we construct the frequency-centered 
light curves giving emphasis in the X-ray and optical bands. 

Results. We show that in cases where 7inax takes low values, the produced X-ray light curves show a plateau as the 
synchrotron component gives its place to the Synhro Self-Compton one in the X-ray band. 

Key words, gamma-rays: theory - acceleration of particles - radiation mechanisms: non-thermal 



1. Introduction 

Gamma-Ray Bursts (GRBs) are attributed to a release of 
a very large amount of energy 10^^ — 10^^ ergs) into a 
small region of space (< 100 km) over a short period of time 
(~ 10 - 10^ s for long GRBs and < 2 s for short GRBs). 
These energetic events have two characteristic radiative sig- 
natures: (i) the prompt 7-ray and (ii) the afterglow emis- 
sion. The detection of high energy photons (e.y ^ 1 MeV) 
implies sources of radiation moving at relativistic speeds 
with Lorentz factors F exceeding 1 00 ([Fenimore et al.iri993l : 
iLithwick & Sarill2001t lPiranlll999h . 



While many issues concerning the prompt emission 
are still open, the afterglow, i.e. the lower energy long 
lasting emission, is believed to arise from the interaction 
of the relativistic ejecta with the ambient matter and 
can be adequately described by the so-called 'standard ' 
model Cr .ccs & Mes zaroslll99l iPaczvnski fc Rhoadsl[T99l 
iMeszaros fc Reeslll997f ). According to this, the relativistic 
blast wave produced from the explosion can energize the 
external medium, i.e. accelerate electrons (and possibly 
protons) to high energies and generate magnetic fields. The 
relativistic electrons radiate by synchrotron and inverse 
Compton radiation which is essentially the observed 
afterglow emission. However, in order to calculate the 
radiated photon spectra, one needs a detailed prescription 
of the electron distribution function and of the magnetic 
field. This is usually done by postulating that the electrons 
have a power law distribution between a minimum (7min) 
and a maximum (7max) cutoff with an overall energy 
density content which is a a fixed fraction (usually denoted 
by £0) of the total post-shock internal energy density 
U while an analogous argument can be made for the 
magnetic field energy which takes a fraction ee of U. A 
significant amount of work has been performed by many 



researchers in calculating the multiwavelength spectra and 
light curves of GRB afterglows either based directly on 
the a b ove prescription (Dermcr & Chiang 1998; Sari ct aL 
1991 iPanaitescu fcMeszarosi 119981: iWiiers fc Galama 



" iPanaitescu fc Kumar 
using different vari- 



'an fc Piran 
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19991; iDermer et al.l I2OOOI: 
200I iGranot fc Sa'rri2002D 01 

ations (iGranot & Kumar' 'm.,- 

Panaites cu et al. 2006; Zhang et al] 20061 : Nousck et al.l 
I2006D . 

In the present paper we focus on the effects that a low 
7max will bring on the multiwavelength spectra and light 
curves of the afterglow emission. This has not been treated 
thus far as it is implicitly assumed that 7max is very large 
and its radiative signature does not contribute to any ob- 
servable band. However, if it has a low value, then a break 
might appear successively in various energy bands of dimin- 
ishing energy as the synchrotron component gives gradually 
its position to the SSC one. This will produce light curves 
which are not any more pure power laws but have more 
complicated shapes. 

The paper is structured as follows. Is §2 we describe the 
principles of the model and discuss, in a qualitative way, 
some of the results. In §3 we quantify the above and we 
derive some analytical relations between the initial param- 
eters, which, when satisfied, will produce different types of 
X-ray light curves. In §4 we make a tentative connection of 
our results to observations. Finally in §5 we conclude and 
give a brief discussion of the basic points of the present 
work. 



2. The Model 

2.1. First principles 

The general framework of the model we present here 
is based on the standard GRB afterglow model, al- 
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beit with some modifications rega rding mai nl y th e 
approac h to the physical problem (jFan et al ] (!2008); 
IPctrop oulou fc Ma stichiadis (200^ - henceforth PM09). 
For the sake of completeness we repeat here its basic 
premises: as the Relativistic Blast Wave (RBW) usually 
associated with GRB afterglows is expanding, it acceler- 
ates by some unspecified mechanism electrons of the cir- 
cumstellar medium to high energies. These are assumed to 
be injected behind the shock front in a volume of radius 
R containing a tangled magnetic field B. The particles suf- 
fer radiative and adiabatic losses, evolving with radius. At 
the same time they emit radiation by synchrotron and in- 
verse Compton radiation. Therefore, at each radius there 
is a coupling between electrons and photons, in the sense 
that the electron distribution function determines the pho- 
ton spectrum and, at the same time, the photons determine 
the electron distribution function through inverse Compton 
losses and, possibly, pair reinjection. The usual procedure 
of approaching the problem is to solve simultaneously two 
coupled kinetic equations for the distribution functions of 
electrons and photons which, when solved, give the afore- 
mentioned quantities as functions of radius and energy. 
The most relevant physical processes which are included 
are: electron synchrotron radiation, synchrotron self ab- 
sorption, inverse Compton scattering (both in the Thomson 
and Klein-Nishina regimes), photon-photon pair production 
and adiabatic losses (for a mo re detailed discussion regard - 
ing the physical p rocesses see iMastichiadis fc Kir^ (jl995l) ; 
iPe'er fc Waxl^ (l200l l. 

To obtain the multiwavelength (hereafter MW) spec- 
trum of GRBs at each radius r of the relativistic blast wave 
one needs to specify the Lorentz factor of the flow r(r), 
the comoving radius of the source R = r/T, the magnetic 
field strength B{r) - determined indirectly through the pa- 
rameter eB, and three parameters related to the electron 
injection, i.e. their total power - determined by the param- 
eter Ce, their slope p of the power law at injection and the 
maximum cutoff of their distribution 7max ~ the minimum 
cutoff 7i„in is defined in terms of the other parameters by 
equation (ITT]) . 

2.2. Multiwavelength spectra 

The approach described above allows one to calculate the 
photon MW spectra self-consistently while it inherently ad- 
dresses questions about whether the electrons are in the fast 
or slow cooling regime. Moreover it can calculate the spec- 
trum in various regimes without resorting in a piecewise 
succession of broken power laws which becomes problem- 
atic in cases where the characteristic frequencies are rather 
close to one another. Finally, it takes into account SSC 
losses which, as PM09 have shown - see also ISari &: EsinI 
()2001|) . can alter significantly the electron spectrum and 
therefore, the photon spectrum, even in the well-studied 
synchrotron regime. 

Some of the above are exemplified in Fig. [T] which de- 
picts a generic case of a MW spectrum in the case of a 
power-law electron injection with 7niax not much greater 
than 7min- As the parameters have been chosen so as the 
radiating electrons to be partly in the uncooled regime, 
the sync hrotron component consists, at least in theory 
(jSari et a l. 1998), of four segments: (1) synchrotron self- 
absorbed part below i^ssa, (2) typical synchrotron low en- 
ergy part for Vssa < v < i^s,min, (3) uncooled part for 
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Fig. 1. Multiwavelength spectrum expected for a case 
where the upper electron cutoff is not much greater than 
the lower one - in the present case 7rnax = 10'* and 
7min,o = 5.6 X 10^. The complete set of the parameters used 
is: = 10^4 ergs, Fo = 400, no = 1 part/cm^, eg = 0.005, 
tc = 0.01 and p — 2.3. Both synchrotron and SSC com- 
ponents of the spectrum appear to be continuously curved 
without any clear power law segments (at least for frequen- 
cies below Vs,m\i\ and I'ssc.min respectively) mainly due to the 
proximity of 7max and 7min- The characteristic frequencies 
indicated in the figure are discussed in the text. 



I's.min < V < Vs,c and (4) cooled part for Vs,c < v < J^s,max- 
One could add here (5) the synchrotron exponential cut- 
off which appears above J^s,max- Between these segments 
power laws of different slopes should appear, however due 
to the proximity of the lower and upper electron cutoffs, the 
turnovers in frequency have been smeared out and the spec- 
trum appears to curve continuously without any clearly de- 
fined power-law regime, at least for frequencies above I's^min- 

The same holds for the SSC component: this is much 
broader, as expected, than the synchrotron one and it 
also shows a continuous curvature. In Fig. [1] three char- 
acteristic frequencies of the SSC component are indicated: 
I'ssc.min = |7minf^s,min, i^sscx = ITc'^s.c and the peak fre- 
quency z/p. Sari & Esin (2001) have shown analytically that 
the peak frequency of the SSC component in the slow cool- 
ing regime is given by I'ssc.c- However, this is not expected 
when a relatively low 7max is taken into account, as in the 
case of Fig. [1] An analytic calculation (which can be found 
in detail in the Appendix) of the SSC peak frequency can 
also be done, in the case of an electron pure power law dis- 
tribution, having limits between 7min and 7max- The SSC 
peak frequency is then given by |7^ax^s,min which equals 
to 1.4 X 10^^ Hz for the example given in Fig. [1] The nu- 
merically calculated peak frequency is however different: 
i/p = 4.4 X 10^^ Hz sa |7c ^'s.min- This simple e xample shows 
that the presence of the Compton logarithm (|Gouldlll979D 
combined with the fact that the electron distribution might 
have at least two breaks, one at 7inin and one at 7c, makes a 
simple analytic approach complicated. Nevertheless, for the 
specific example one can clearly see in vF^ units a rising 
part, a broad peak and a declining part. 
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2.3. Light curves 

Based on the above it is straightforward to construct hght 
curves at various frequencies. However, it is worth describ- 
ing quahtatively a case where cutoff effects in the syn- 
chrotron component appear due to a low value of 7max- 
Assume a certain observing window between two frequen- 
cies. Assume also that the initial parameters are such as 
initially the flux in this window to be dominated by the 
synchrotron component. As time evolves in the observer's 
frame, he/she will first observe the various parts of the syn- 
chrotron c omponent passin g through it, as first was sug- 
gested by ISari et all (|1998[ ). Therefore at some point in 
time, which we will call fbr.i, the combination of F, B and 
7inax will be such that the flux in the observing window will 
be dominated by the exponential cutoff of the synchrotron 
component. This will result in a natural steepening of the 
light curve. In the hypothetical case where the SSC com- 
ponent is absent, the observer would have seen an ever in- 
creasing steepening of both the spectral index and of the 
light curve until the fiux would drop to very low levels, be- 
low the sensitivity of any instrument. In reality, however, at 
some level the SSC component must appear in the observ- 
ing window and start dominating the fiux there. This would 
naturaUy result in a fiattening of the light curve as the ob- 
server starts sampling photons not from the synchrotron 
but from the SSC component. Depending on the relation 
between jmin and 7max, at time t = tbr,i the SSC compo- 
nent could be, broadly speaking, in any of its three spectral 
regimes (i.e. rising, peak or declining). As we elaborate in 
the next section, if the SSC component becomes dominant 
while still early in its rising part, the light curve will show 
a sharp turnover which will result in a shallow decline of 
the flux, i.e a plateau. If, on the other hand, the transition 
from the synchrotron to SSC occurs while the latter is close 
to the peak or during its declining part, the light curve will 
show a more gentle flattening. 

Figures [2] and [3] present the results of a run showing an 
X-ray plateau. Figure [2] shows snapshots of the MW spectra 
at five different times, obtained for a GRB at z = 1 with 
Eo = 10^4 ergs, Fq = 400, uq = 1 part/cm^, eg = 0.001, 
eo = 0.025, p = 2.3 and 7max = 2.24 x 10''. One can see 
that at early times the X-ray band is dominated by the tail 
of the synchrotron component which in about t ~ 10"^ s 
has given its position entirely to the SSC one. This trend 
is repeated once more at much later times (i ^ 5 x 10^ 
s) in the optical band. Figure [3] shows the corresponding 
optical and X-ray light curves. The latter shows a clear 
plateau which gradually steepens due to the concavity of 
the SSC component. On the other hand, the optical light 
curve shows the 'standard' unbroken power-law behavior 
until t ~ 5 X 10^ s and flattens off slightly at even later 
times, as the SSC component takes over. 

The above results show that in order the X-ray light 
curves to show plateaus, at least within the context of our 
model, the SSC component should be rather flat in the X- 
ray energy range at the time the synchrotron component 
decays due to its exponential cutoff. Since the above condi- 
tions implicate the magnetic field B{r), and the lower and 
upper cutoffs of the electron injection, 7min and 7max re- 
spectively, one can quantify the above conditions using the 
parameters of the standard afterglow model. We proceed 
to do this in the next section. 




Fig. 2. Multiwavelength spectra at observer times from 10^ 
s to 10^ s (top to bottom). For the parameters used see text. 
X-ray and optical windows that correspond to the observing 
energy ranges of XRT (0.3-10) keV and UVOT (170-650) 
nm respectively, are also shown. 
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Fig. 3. X-ray (solid line) and optical (dashed-dotted line) 
light curve corresponding to the same case as in Fig|21 The 
contribution of the synchrotron (dashed line) and SSC (dot- 
ted line) components to the total X-ray flux are also shown. 
The optical light curve is shifted by -1-2.5 units in logarithm 
for reasons of better display. 



3. Effects of 7max 

The upper limit of the electron distribution 7niax has not 
been so far considered in GRB afterglow models as a dy- 
namic parameter. In this section we will show how a rela- 
tively low 7i„ax can affect the MW spectra and the corre- 
sponding X-ray and optical light curves. We assume that 
the afterglow is produced by an adiabatic relativistic blast 
wave decelerating while interacting with the ISM. During 
this phase the evoluti on of the RBW is describe d by the 
self-simflar solution of iBlandford fc McKe^ ()1976[ ): 



-3/2 



ff r > Rd, 



(1) 
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where 
Rd — 



1/3 



Using the relation 



(2) 



(3) 



and neglecting the term = compared to /^^ 2cr'^{r') 
we find that 

r « (8ctrlRiy/\ 



(4) 



For the strength of the magnetic field (in the fluid frame), 
we adopt the usual form 



B{r) = y^32TTnonipeBcT{r). 

The observed synchrotron frequency corresponding to "fn 
is given by: 



(5) 



r- 



eB 
27rmp 



-7' 

/ max 



1/4 



/max" 



^3/4 



(6) 



We consider an X-ray window that corresponds to the ob- 
serving energy range of XRT, i.e (0.3 — 10) keV. We are 
interested in the time t^^^ ^ when photons belonging to the 
synchrotron exponential cutoff which forms above t's,max 
(see Fig. 1) cross a characteristic frequency of the X-ray 
band, say I'x ~ 6.3 x lO^*" Hz. So we set Vs.coS = ^^^s.max 
where A is a numerical factor of order 5—10 which deter- 
mines how deep into the exponential cutoff the particular 

(x) 

synchrotron photons are. Then the expression for t^Jf^ ^ be- 
comes 



t. 



(x) 
coff ,s 



V 327r3c5mf 



1/3 



-4/3^4/3 8/3 . 
X /max 



(7) 



Demanding that the characteristic time i^^jj , 
interval of the general form: 



tk < t 



< tk+l 



lies in a time 



(8) 



where tk — 10*^ s, equations ([T]) and (0 
first constraining relation for 7max: 

gi{Eo,no,£B) < 7max < g2{Eo,no,£B), 



I combine to give the 



(9) 



where .gi,2(i5o, no, ee) = Ci,2 10^'=/' A-^^^E^^i^n^'^/'e^'ll 
with Ci = 4.1 X 10^ and C2 = 9.7 x 10^. One should 
keep in mind that relation ([8]) is valid only for values of 
the real variable k which ensure that time tk is larger 
than the deceleration time td- Here and throughout the 
text, the convention Qx = Q/10^ has been adopted in 
cgs units. It is worth noting that for typical values of 
^'0,54 = no,o = eB,-2 = 1 and k = 2 the maximum electron 
Lorentz factor lies between 8.1 x 10"^ and 1.9 x 10'', which 
is relatively low. Although the above double inequality in- 
volves three free parameters of the model, the dependence 
on two of them, i.e on Eq and no is very weak. 



The SSC component peaks at a characteristic frequency 
which depends on the Compton logarithm first introduced 
by Gould (1979). In the simple case where the electron dis- 
tribution is given by a pure power law between 7niin and 
7max and the scatterings happen in the Thomson regime, it 
is found that the peak frequency of the SSC component is 
^'p = |7max^s,min (scc Appendix for a detailed calculation). 
In a generic case where the electron distribution shows a 
cooling break, the calculation of the SSC peak frequency 
is more complicated (see discussion in §2.2). However, in 
all cases the minimum frequency of the main SSC branch, 
1-6 ^^ssc,min — iTmin'^s.min is bclow (slow cooliug regime) or 
at least equal (fast cooling regime) to the peak frequency. 
For that reason we choose to quantify the lower energy part 
of the SSC component as 6i'ssc,min, where 6 is a numerical 
factor of order 0.001 — 0.01. In order to proceed we need to 
use an expression for 7niin- So far in GRB afterglow mod- 
els, 7max was not treated as a 'dynamic' parameter in the 
sense, that its signature would not be observed in the X-ray 
energy range and below, as it was taken to be much larger 
than 7niin- For that re ason the approximate expression for 
7min (|Sari et al.lll"998h 



approx nip p ^T-i/ \ 



' rric p — 1 



(10) 



was safely used. However, in our work where we also exam- 
ine cases with 7inax only a few times greater than 7,nin, in 
the numerical code we use the accurate expression which is 
solution to the equation: 



2—p 2—p 

7m in Tmax r II!lE.Vfr\ 

' p 1~P ~ "■'P m ^ 
— 7max '"'(^ 



7, 



(11) 



where fp 



As the solution to the equation above 



p-i- 

has no explicit analytical form, we proceed first to find a 
constraint of 7max using the approximate analytical expres- 
sion given by eq. ((TU]). For the low energy part of the SSC 
component we can now write: 



6m 



SSC, mm 



2eb 



■^327rnoTOpeB ( e-c — fp ) 



/ 3£;o 



\20487rnompC" 



3/4 



-9/4 



(12) 



In order to have a flattening of the light curve after the 
first break, this low energy part of the SSC component 
must appear in the X-ray band about the time when the 



synchrotron one decays. In complete analogy to t 

observed time t^l at which the frequency hv, 
the X-ray band is determined by the relation 



(x) 



the 



SSC, mm 



coff ,s' 

enters 



f(x) 



4/9 f _ 



2eb 



\37rmg 
3£;o 



4/9 



20487rnompC 



(327rnoTOpeB)^^^ 

1/3 



rup 

Cc Jp 



16/9 



(13) 
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The ratio of the two characteristic times is given by: 



t'SSC 
''C0ff,E 



"6^ 



4/9 



^4/9^-4/3 -8/3 
'max 



noCB I 



2 X 1016 6*/^A-4/3^^8/3 ( j/^^e 



4/9 



noCB 



4/9 

(14) 
(15) 



fx) 

When the two timescales i^^jj ^ and t, 



ssc are of the same 
order, a break will appear in the X-ray light curve, as the 
synchrotron component gives its place to the low energy 
part of the SSC one. Allowing 



0.2 < 



t 



coff.s 



< 1.0 



(16) 



we obtain one more constrain for 7max, i-e. 

hlino, tc, CB, /p) < 7max < ^2(n-o, Co, eb, /p), 

where the functions hi 2 are defined as /11.2 



(17) 



2/3 2/3 



-l"'0,0 



^^^Cr^/j with = 6 X 10^ and 



iir2 = 1.1 X 10^. The above relations make the tacit assump- 
tion that the scatterings occur in the Thomson regime, at 
least for the lower energy part of the SSC spectrum. For 
that reason, we impose one more constraint on 7max de- 
manding that at least for t ~ ^ the assumption of scat- 
tering in the Thomson regime is valid. For that we use the 
dimensionless factor 



7n 



he ( 3i?oeB 



1/2 



-3/2 



(18) 



where primed quantities are measured in the comoving 
'•coff s ^ parameter is given by: 



. Att = 


y:(x) 
''coff.s 


hf^ml 


/327r 







1/2 



4 4^ ^7mL("oeB) 



-1/2 



Thus, when a; < 1 we find that 
> 50 A-^/^n, 



7. 



1/8.3/4 3/4 -1/8 
0,0 Jp *'o,-l*=B,-2- 



(19) 



(20) 



The above relation for typical parameter values poses only 
a weak constraint on 7max- The same holds even if the ac- 
curate expression for 7min was used. For that reason from 
here on we will not take into consideration the constraining 
relation (QUI) . 

Having determined the constraining relations for 7„iax 
using the approximate form of 7min, we can now estimate 
the corrections introduced, after taking into account the 
accurate expression of it. For the purposes of our analytic 
analysis we model it as: 



7n 



approx 
mod7min ' 



(21) 




Fig. 4. Log-Log plot of /mod as a function of radius r (in 
units of the initial radius Rq = lO^** cm) for two values 
— : 1.32 (solid line) and 33 (dashed line). 



of the ratio 



Other parameters used are: Eq 
no = 1 part/cm^, €b = 0.001, eo 
deceleration radius is also shown. 



10^-* ergs. To = 400, 
0.01 and p = 2.3. The 



reasons, against radius r for two values of "^'""'"^ in Fig. 01 It 

is evident that /mod —J' 1 in the limit of 7max >> 7min,o, as 
expected. As our analysis holds for the deceleration phase 
of the RBW, we can simplify the calculations further by 
averaging the function /mod over f = logr for r > Rd'- 



mod 



7n 



7min,0 



r>Rd 



drf mod ( f 



Ir 



r>Rd 



df 



(22) 



This average modifying parameter can be used in eq. (1211) 
instead of /mod- Thus, constraining relation (1171) becomes: 



hi < 7max < h2, 



(23) 



where ft.1,2 ~ /mo'cf^i.s- Table [T] shows some indicative 
values of the correction introduced in relation for 
Ea = 10'^^ ergs, Fq = 400, no = 1 part/cm^, ee = 0.001 
£0 = 0.01 and p = 2.3. 



log 7max 


log (7max/7min,o) 


1 log /mod 


4.15 


0.495 


0.087 


4.35 


0.787 


0.072 


5.35 


2.00 


0.030 



Table 1. Indicative values of the correction introduced in 
the constraining relation (P^ . For the parameters used see 
text. 



Using eq. ^ and ([^ we can plot 7max versus £0 for 
fixed values of Eg, no, p and eb, creating a parameter space 
shown in Fig.[5j The other parameters used are the same as 
those used in Table [T] above. The curves defined by eq. (O 
and (j23p create distinctive areas on the parameter space. 
For initial values chosen from the striped area, eq. dTTI) 
has as solution 7min,o = 7max with no physical meaning. In 
that sense the striped area is not permitted. The horizontal 
zones labeled by A: = 2 and k — 3 are related to cases 

(x) 



where the term /mod is a function of the radius r and of where the transition time t^^jj ^ lies between 10^ — 10^ s and 
the ratio ; 7min,o is the initial minimum Lorentz factor 



of the electrons. Function /mod is plotted, for illustrative 



coff ,s 

lO'^ — 10^ s respectively. A choice of fc < 2 would correspond 
to breaks occuring at even earlier times. We call a 'plateau 
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strip' the tilted zone, as a choice of a pair (eoj7max) from 
it leads to X-ray light curves which show a shallow decline 
phase, i.e a 'plateau'. Moreover, if this pair is chosen from 
the intersection of the 'plateau strip' with a horizontal zone 
labeled by /c = 2 for example, the plateau phase will begin 
at a time between 10^ — 10'^ s and so on. Finally, a choice 
of pairs from the area outside the tilted zone leads to X-ray 
light curves with a single change of slope. 

For a better inspection of the parameter space, we have 
chosen pairs from four distinctive regions of the diagram 
and show the corresponding X-ray and optical light curves 
in Figs. El and [5] As we move from point 1 to 2, i.e by 
increasing while keeping 7niax constant, the break times 
(see Fig. [5]) remain constant while the light curve shapes 
exhibit a transition as the decay becomes flatter leading to 
a plateau phase. From point 2 to 3, the break is shifted 
towards later times while the light curve becomes, once 
again steep. This is to be as expected since point 3 lies 
outside the 'plateau strip'. Finally, as we move from point 
3 to 4 the slope of the light curve after the break becomes 
flat and a plateau is produced. This behavior is the same 
with the already discussed transition from point 1 to 2, with 
the sole exception that now the break occurs at later times. 
Note that for all the numerical runs presented in the present 
work, we use the accurate expression for 7min by solving 
eq. pT|) . For that we have included in our numerical code a 
subroutine that utilizes a c ombination of the bisection and 
Newton-Raphson Method (jPress et al.lll992D . 




Fig. 5. Initial value parameter space of 7max and eg for 
Eq = lO^"* ergs, Fq = 400, no = 1 part/cm3,eB = 0.001 
and p = 2.3. The numerical factors chosen here are A = 
8,b — 0.001 (see text for their definition). Dotted lines 
correspond to the constraining relation solid lines to 
and the dashed line sets the boundary of the striped 
area, which is not permitted as it it leads to 7max = 7inin,o- 
The horizontal zones labeled by fc = 2 and k — 3 cor- 
respond to transition times t^^g ^ lying in the time inter- 
vals 10^ — 10"^ s and 10^ — lO^'s respectively. The 'coordi- 
nates' (eo,7max) of the points 1 to 4 marked on the plot 
are: {(0.005,1.41 x 10*), (0.01,1.41 x lO'*), (0.01,2.24 x 
10^), (0.025,2.24 X 10'')} respectively . The corresponding 
X-ray and optical light curves of the aforementioned points 
are presented in the following Figs. [5] and [SI 




Fig. 6. X-ray light curves for different set of parameters 
corresponding to the points from 1 to 4 (bottom to top) 
marked on panel (b) of Fig. [S] The first break of the X- 
ray light curves 1 and 2 occurs at around 125 s, while for 
the light curves 3 and 4 occurs later at i ^ 650 s. For 
clarity reasons light curves 1, 3 and 4 are plotted with an 
offset of —0.2, -1-0.4 and -1-0.5 in logarithmic units of flux 
respectively. 



One could also make some comments on the effects of 
the other free parameters, which were assumed fixed in 
Fig. [5l will have on the parameter space. For example, a 
possible change of Eq will affect only the dotted lines in 
Fig. [SJ The effect will not be of great importance as the in- 
equality ([9]) has a very weak dependance on Eq. In general, 
an increase of ee will shift the 'plateau strip' downwards. 
Although the parameter Fq does not appear explicitly in 
the constraining relations (HJ and ((23|) has an effect on the 
relative positions of the striped area and the 'plateau strip'. 
Finally, the external number density uq may not affect 
severely the appearance of the parameter space, but has an 
important effect on the ratio of the synchrotron to the SSC 
flux. We note that as it was shown in PM09, an increase 
of Uq makes the afterglow more Com pton dominated. This 
possibility has also been discussed bv lPanaitescu fc Kumail 
(2000) and Sari & Esin (2001). A denser external medium, 
increases in vF^, units the X-ray flux over the optical. Such 
high values of external density are required by our model 
in cases where the X-ray light curve shows a shallow phase 
and the ratio Lx/iopt > 1- Figure [7] shows such a case. The 
parameters used are Eq = 10^"^ ergs, Fq — 100, uq = lO'^, 
ee = 10-5^ = 0.032 and 7max = 8.2 x lO''. For t > 10^ 
s, when the spectral eivolution in the X-ray energy band is 
not very important, the optical flux lies approximately one 
order of magnitude below the X-ray one. 

As the RBW decelerates the synchrotron and SSC emis- 
sion from the forward shock becomes weaker and softer (see 
Fig. [2]). As a result, the synchrotron cutoff will enter at 

some instant i|;°g*s in the optical band. It is straightfor- 
ward to show, under the assumptions of the problem, i.e. 
when B oc r~'^/^ and 7max independent of r, that this is 



related to t 



.(opt) 

'•coff.s 



(x) 
coff ,s 



''opt 



by: 



4/3 



'•coff.s 



(24) 
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Fig. 7. Multiwavelength spectra at observer times 10'^ s 
and 10^ s (top) and light curves for the X-ray and optical 
band (bottom) for an SSC dominated afterglow. For the 
parameters used see text. The X-ray and optical windows 
are as given in Fig. [21 

or, equivalently, by 

t^l « 2 X 10^ ^(^1,,,, (25) 

where a typical optical frequency I'opt = 2 x 10^^ Hz was 
used. For the cases represented by points 1 and 2 in Fig. [5] 

one finds t^^g ^ « 125 s (see also Fig. 15]). Thus, a break in 

the optical light curve is expected at around 2.5 x lO'^ s, 
which is confirmed numerically and it is shown in Fig. |51 
The other two optical light curves in the same figure show 
a break at correspondingly later times (> 6.3 x 10^ s). We 
should also note that the break of the optical light curve 
will not neccessarily be of a 'plateau' type. This is due 
to the fact that the emerging SSC component has become 
steeper since the time it has emerged in the X-ray window 
(see Fig. 2 for the time evolution of the SSC component). 
A natural outcome of our model is that the optical light 
curves do not exhibit a break during the plateau phase 
of the correspon ding one in X-rays, as it was first noted 
by iFan fc PiranI (2006, ) . Evidence for chr omatic breaks in 
X-ray light curves were also discussed in iPanaitescu et aD 
((200 6'). 

Chromatic breaks in the optical and X-ray light curves 
can be produced by our model with an expected time dif- 
ference given by eq. (^5]) . As we discuss in the last section, 
this is a strong constrain which reflects the validity of the 
model's assumptions. 

4. Basic results 

It is well known that the X-ray and optical light curves of 
GRB afterglows sho w a wide range of behaviors. Especially 
in the X-ray regime (|Nousek et al.ll2006l: IEtous et al.ll2009l ) 
the behavior can be quite perplexing (see Fig. ^ with many 
cases showing a plateau which is quite difficult to interpret 
within the context of the standard model. The analysis per- 
formed in §3 can tentatively be of some relevance to these 
observations. Using, for instance the example set by Fig. [5] 



-8 I 




2 3 4 5 5 

log t (S) 

Fig. 8. Optical light curves for different set of parameters 
corresponding to the points 1 to 4 (bottom to top) marked 
on panel (b) of Fig. [5l For clarity reasons light curves 1, 3 
and 4 are plotted with an offset of —0.4, +0.2 and +0.5 in 
logarithmic units of flux respectively. 

as a guide, we can say that plateaus appear as long as 
the values of 7max and te are chosen from the region de- 
fined by the two tilted lines. Moving inside this region from 
the lower left to the upper right, plateaus appear at pro- 
gressively longer times. Furthermore, a choice of the initial 
parameters outside this region (for example, from the left 
top corner of Fig. [SJ leads to afterglows without a plateau 
phase. Fig. [TU] shows different types of X-ray light curves 
obtained using our numerical code corresponding to points 
from different regions of the parameter space of Fig. [5] Light 
curves in panels (a) and (b) correspond to points (4) and 
(1) already shown in Fig.[5l while light curves of panels (c) 
and (d) are obtained using (ec, 7max) = (0.0032,6.3 x 10^) 
and (0.01, 10^) respectively^ A tentative comparison of our 
model light curves to those of Fig. [9] can be made. 



V a) 









Time 



Fig, 9. Schematic diagram from lEvans et all (|2009D show- 
ing the different light curve morphologies observed. 

As the model presented here produces multiwavelength 
spectra at each instant, we can use it to calculate the evo- 

^ All other parameters used are the same as in Fig. [S] except 
for the case presented in panel (c) where To = 800 was used. 
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Fig. 10. Different light curve morphologies obtained using 
our numerical code corresponding to points from different 
regions of the parameter space shown in Fig. [5] Light curves 
from each panel can be tentatively compared to the corre- 
sponding ones of Fig. [3] For the parameters used see text. 
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Fig. 11. Time evolution of the hardness ratio for each of the 
cases presented in Fig. [TO] For the definition of the hard- 
ness ratio used see text. The insert in each panel shows the 
evolution of the corresponding photon index. The shaded 
areas in the inserts imply that during this period our model 
spectra cannot be approximated by a single power law and 
therefore a specific photon index could not be attributed to 
the spectrum. 



lution of the expected X-ray hardness ratio defined as the 
ratio of count s in the 1.5- 10 keV to the counts in the 0.1- 
1.5 keV band (|Evans et ah 2009. 20I0i)- Figure [H] shows the 
time evolution of the hardness ratio for each of the example 
cases shown in Fig. [TUl Time evolution of the correspond- 



ing photon index is also shown in the inserts of Fig. Illl 
whenever the spectral shape allows its viable calculation at 
the particular time (for a more detailed discussion on the 
shape of our X-ray model spectra see §2). 

For X-ray light curves with a distinctive 'plateau' phase, 
we find that the spectral evolution shows a characteristic 
trend as the X-ray window is first dominated by the syn- 
chrotron and later by the SSC component. This can be seen 
in panel (a) of Fig. [TT] At very early times both the soft and 
hard X-ray bands are dominated by the synchrotron pho- 
tons - however the hard band is affected first by the syn- 
chrotron cutoff and this has as a result the decrease of the 
hardness ratio. During this phase the spectrum in X-rays 
is shaped by an exponential cutoff (synchrotron emission) 
and an emerging flat power law component (SSC emission). 
Thus, it cannot be simply approximated by a single power 
law and 'labeled' by a photon index (shaded area in the in- 
sert). At later times, the SSC component starts appearing 
in the hard band while the decreasing synchrotron compo- 
nent dominates the soft one, resulting in an increase of the 
hardness ratio. Finally, at even later times both bands are 
dominated by the SSC component, whose low energy part 
can be approximated by a flat power law, and due to its 
gradual steepening the hardness ratio appears to decrease 
gently. 

In cases where the X-ray flux decays as a power law with 
time, as in panel (d) of Fig. llOi we find no significant spec- 
tral evolution. The photon index is approximately constant 
almost for three or four decades in time, as the power law 
segment of the synchrotron component dominates until late 
times in the X-rays (see insert in panel (d) of Fig. [TT|) . The 
other two cases presented in panels (b) and (c) of Fig. [TT] 
lie somewhere inbetween the two aforrementioned example 
cases. 

Although the qualitative evolution of the photon index 
with time is a robust feature of our model , its specific value 
depends on the value of the other model parameters, such 
as the slope of the electron energy spectrum. In all our runs 
we have used a typical value of p = 2.3. 

5. Summary/Discussion 

In the present paper we have investigated the role that the 
upper cutoff of the electron injection can play in the evo- 
lution of the multiwavelength spectra and light curves of 
GRB afterglows. For this we have solved self-consistently 
the kinetic equations that govern the electron evolution 
and photon rad iation as a function of distance (see also 
iFan et al.l ()2008[ ) and PM09). This approach can address 
successfully the effects of the electron cutoff radiation on 
the light curves. 

We have shown that depending on the adopted value of 
7max the X-ray lightcurves can show one of the following 
behaviours: 

1. In cases where 7max is not much greater than the lower 
cutoff 7min, the X-ray light curves show three distinct 
phases. First a fast drop phase which corresponds to 
the exponential cutoff of the synchrotron component. 
Then a plateau phase which is caused from the gradual 
dominance of the SSC component over the decaying 
synchrotron. Finally, a more gradual power-law decay 
which corresponds to the normal evolution of the SSC 
component. The analytical approach used in §3 and 
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summarized in Fig. [S] shows that, depending on the 
initial parameters, there might be a narrow strip in 
phase space that produces X-ray plateaus. In all cases 
we found that 7max should be no more than a factor of 
10 larger than 7niin- 

2. If 7max does not satisfy the above condition but it 
is still close to the 'plateau strip', then the X-ray 
afterglows do not show a plateau but simply a change 
in the power law decay, i.e. the three phases degenerate 
into two. 

3. In cases where 7max ^ 7min we find that the X-ray 
afterglows will be dominated until very late times by 
the synchrotron component, i.e. we obtain the standard 
afterglow picture. 

These trends are exemplified in Figures [S] and [TUl As 
far as the fiux decay in the optical regime is concerned, 
the present model predicts that the optical light curve will 
mimic the X-ray one with two major differences (see Fig. [5]). 
First, the break in the optical light curve will come at much 
later times (see ea. ips]) ) and second, after the break, the op- 
tical light curve will not neccessarily show a plateau. Note 
however, that both of these statements are based on the as- 
sumption that 7inax remains constant throughout the RBW 
evolution. 

In this respect, 7max emerges as one of the important 
parameters of the afterglow evolution as its choice can con- 
trol critically the behaviour of the X-ray light curves. 

The evolution of the X-ray hardness ratio and spectral 
indices were presented in Fig[TTJ Our model derived hard- 
ness ratio shows a characteristic signature which is com- 
patible with observ ations at least du ring t he early stages 
i|Butler fc Kocevskil[2007t iLiang et al.ll2007D . On the other 
hand it is still inconclusive regarding the late stages (P. 
O'Brien - private communication), as the transition from 
the plateau to the normal decay shows an evolution in the 
hardness ratio by a factor of two, which is not the cas e 
for several i ndivid ual GRBs (e.g. IVaughan et al.l ()2006D : 
iLiang et all (|2007D l However, as preliminary calculations 
have shown (Petropoulou, Mastichiadis & Piran (2011) - 
to appear in the proceedings of the 25*^^ Texas Symposium 
held last December in Heidelberg) , this is greatly relaxed in 
the case where 7max is allowed to increase with radius. As 
the study of such cases is beyond the scope of the present 
paper, we will treat this in a forthcoming publication. 

It is important to emphasize at this point that the close 
relation we found between X-ray light curves exhibiting a 
plateau phase and a 7max not much greater than 7inin, is 
what actually differentiates our model from other works 
(e.g. long lasting energy injection into the forward shock 
- refreshe d shock models (Zhang et al. 2006), late prompt 
emission (iGhisellini et al.i ,2007: .Ghisellinl .2008il. geomet - 
rical effects (|Eichler fc GranotI 120061: iGranot et al l [2006[l . 
non-standard deceleration of the bulk Lorentz factor due 
to the Compton- drag force ([Mas tichiadis & Kazanas 2009;) , 
dust scattering (jShao et al.l l2008l ) and others). Some re- 
cent PIC simulations show that most of the particles are 
accelerated into a relativistic Maxwellian, while a small 
fraction of them is injected to a power law high energy 
tail, whose high energy cutoff is approximate ly only one 
order of magnitude larger than the low one ([Spitkovskvl 
l2008l:ISironi' fc SDitkovskvl[2009l) . However, this is not a fi- 



nal result, as the present numerical simulations have not 
reached yet a steady state, where 7inax is expected to be 
larger than the one found so far. It is interesting to note 
also that when we set 7max not much larger than 7i„in the 
qualitative behavior of our results is similar to the one ob- 
tained when one replac es the power law injection w ith a 
relativistic Maxwellian (jGiannios fc Spitkovsk^l2009h . Ah 
other parameters used in our work are the same as in the 
standard GRB afterglow model. 

Spectral evolution during the steep early phase, which 
is observed i n a s ignificant number of GRB afterglows (see 
I Zhang et al.l (2007) for a systematic study of 44 steeply de- 
caying X-ray afterglows), is an inevitable outcome of our 
model, as in our present work the early steep decay of the 
X-ray light curve is atrributed to the emission from the 
external shock. Moreover there is some observational evi- 
dence for a possible smooth co nnection of the early a fter- 
glow to the prompt emission (jBarthelmv et al.ll2005D . In 
such a case the present wor k could be seen as an extension 
of the supercritical model (|Mastichiadis fc Kazanasll2009f) 
to the afterglow regime. Alternatively it can result from 
a gradual transition towards the end of the prompt phase 
to an external shock emission. Evidence for such a tran- 
sition is seen in several examples where the extrapolated 
BAT light curve is not connected to early-time XRT light 
curve (Tagliaferri et al. 2005). This fact combined with a 
strong spectral evolution at early times suggests that the 
two emissions are produced by either different mechanisms 
or in different regions. Such a transition from the prompt to 
the afterglow emission could be explained wi thin the con- 
text of the int ernal /external shock scenario (iPiranI (Il999l 
200^- see also iDerme 3 (|200l for a discussion on the pos- 
sible scenarios). 

Finally we would like to present potentially problematic 
points of our model - some of these could be addressed with 
future observations and further analysis. 

1 . In the case of an afterglow with a plateau phase we find 
that the X-ray spectrum cannot be represented by a 
simple power-law during the steep decay-early plateau 
phase. The spectral shape at this stage can more 
accurately be described by a steep component coming 
from the synchrotron cutoff plus a flat power-law 
coming from the emerging SSC component - see Fig. [3] 
at early times. 

2. Some preliminary efforts in fitting lightcurves and 
hardness ratios □ have shown that, in some cases, 
our model can successfully reproduce both (e.g. GRB 
060512 - see Petropoulou et al. 2011). However, in other 
cases the model can fit successfully only the lightcurve 
while the HR fit is poorer (e.g. GRB 050713B). 

3. If the light curve shows an abrupt break (either a 
plateau phase or a single change in the slope decay) 
and Lx/Lopt > 1, then our model requires high values 
of the external density {n > 100 cm^'^). This could be 
problematic in cases of GRB afterglows for which low 
values of Nh are derived (Schadv ct al. 2007). We note 



^ Although a successfully fit of the HR could mot neccessarily 
mean a succesful spectral fit, we preferred to model the former 
than the latter as in many cases our X-ray spectra cannot be 
reproduced well by a single power-law - see also the previous 
point. 
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however that, in all other cases the above constrain of 
our model is relaxed. 



the Compton synchrotron logarithm we can safely proceed 
using the spatial averaged photon density 



4. For those GRBs where the early steep decay is smoothly 
connected to the prompt emission phase, the model re- 
quires that the late prompt emission is already domi- 
nated by an external shock emission. This arises natu- 
rally in the external shock model or it requires a transi- 
tion, during the late prompt phase , from in the internal 
to external shock. 

Concluding we could say that the consideration of the 
electron distribution's upper energy cutoff as another free 
parameter in the standard afterglow model brings many 
interesting features in the light curve / spectral behavior 
of GRB afterglows. All these are related to the evolution 
of the relativistic electron distribution and eventually point 
out at the acceleration mechanism at work. We plan to deal 
with this issue in a forthcoming publication. 
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ns{e) =<hs> t 



(A.4) 



where < rig >= J dr r'^hs{r). We assume further that 
the Inverse Compton emissivity is given by a i5 - function 
centered at the mean energy of an upscattered synchrotron 
photon of energy e 



jie(ei;7,e) = ^ei'5(ei-4/37M 



(A.5) 



where A is a normalization factor. Then the total SSC 
power per unit energy emitted is found by 



-'^ic(ei) = y * "-s(e, J ^7 A^c(7)iic(ei;7, e) 



(A.6) 



In the above equation iVo — K^j^^ is the total number 
of electrons per Lorentz factor 7 in the source. The nor- 
malization factor Kc is related to fee through the integral 
Att J dr r'^ko(r). In order to simplify further the calculation 
of the integral in (jA.6|) we use the average photon density 



Jic(ei) = AK^ <ns> 



dee 



(A.7) 



where emin, Cmax are the effective minimum and maximum 
energies of the synchrotron photons and 



Appendix A: Synchrotron Self-Compton spectrum 
of a power-law electron distribution 

We assume a spherically symmetric source of radius R with 
a magnetic field B randomly oriented. Synchrotron radia- 
tion is produced by an isotropic electron distribution, of the 
form 



fcc7 



7n 



< 7 < 7n 



(A.l) 



where p is the exponent of the power-law spectrum, 
7min,7max are the cutoffs of the distribution and is a 
normalization factor that determines the electron density 
and depends on radius r in the case of an inhomogeneous 
source. The synchrotron photon production rate per unit 
energy and per unit volume is then given by 



2e^ 



f 3e 



dN. 



dedt TTich'^c'^ KAnmcC 



a{p) kcB~ e~ 



(A.2) 



w here a{p) is a combina t ion of F - functions (see eq. (4.60) 
of 'Blu menthal fc Gouldl ()1970D ). Equation (|X2|) holds for 
photon energies not close to the low- and high- energy ends 
of the spectrum - es7min ^^'^ £s7max respectively, where 
eg = c ■ "^^^ number density of synchrotron photons in 
the source is given by: 



ns{e,r) = tn 



dNs 
dedt 



ns{r)e~ 



(A.3) 



where tj^ is the crossing time of the source and fig 
ergy independent factor on the right hand side of eq. 
The dependence on radius comes from the quantities fee and 
B. The expl icit functiona l form of the photon density can 
be found in iGouldl (|l979l) . However, for the calculation of 



dj 7 PeiS{ei - 27^e) 




(A., 



for 37inin^ < ci < -^Jmax^- Strictly speaking, the integral 
equals to zero for any other value of ei . However, if one uses 
the complete expression of the Inverse Compton emissivity, 
finds that outside this energy range the intensity is highly 
reduced but not actually zero. The lower limit of integration 
in (IA.8I) is determined by the kinematics of the compton 
scattering and the lower cutoff of the electron distribution, 
i.e 



7min = max[7i„in, (3ei/4e)^/^]. 

The upper limit of integration is given by 

7max = min[7,nax,nioC^/e] 



(A.9) 



(A.IO) 



which takes into account the effect of the Klein-Nishina 
cutoff. For reasons of simplicity we proceed to the calcula- 



tion of the SSC spectrum assuming that % 



7n 



and 



7n 



7n 



Appendix B: The Compton-synchrotron logarithm 

The integration over the synchrotron photon distribution 
(see eq. (IA.7I) . (|A.8p ) leads to a factor 



InS = In 



(B.l) 



called the Compton synchrotron logarithm (|Gouldl[l979[ ). 
This quantity takes into account the effective minimum 
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and maximum energies of synchrotron photons, which con- 
tribute to the principal branch of the SSC spectrum, i.e to 
upscattered photons with energies between ~ es7min ^^'^ 
~ is^tinx- The actual extent of the synchrotron power law 
segment is 



es7min < e < es7max 



(B.2) 



where = ^ . In the previous section we showed that 
only for energies of the upscattered synchrotron photons 
between 



(B.3) 



the SSC spectrum differs significantly from zero. From eq. 
(|B.2I) . (|B.3p we find e < es7max and e < 3ei/47^i„. 



(B.4) 
(B.5) 



1. If 3ei/472jji„ < EsT^ax then one sets 
emax = 3ei/47^i„ 

Cmin — £s7min- 

Thus, 

V - 3£i 1 4 . . 1 2 2 

4e.7mi„ 3 3 

2. If 3ei/472^i„ > e,7^,, then 



(B.6) 



(B.7) 



However, synchrotron photons with energy es7^in can- 
not be upscattered to energies ei > |es7min7max- 1'^ this 
case the effective minimum energy of the synchrotron 
photons is given by 



/max 

Thus we find 

J. ^ 4e,74 4 ^ _ 2 „2 



3ei 

Summarizing, 

(p-3) 

eiJicae, ^ 



2^s7min7max < ^1 < 2^s7inax 



(B. 



(B.9) 



In 



In 



(45") 



< ei < 3es7mi„7 



mm /max 
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